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Typical radio pulsars are magnetized neutron stars that are born rapidly rotating and slow down as 
they age on time scales of 10 to 100 million years. However, millisecond radio pulsars spin very rapidly 
even though many are billions of years old. 1 The most compelling explanation is that they have been 
"spun up" by the transfer of angular momentum during accretion of material from a companion star 
in so-called low-mass X-ray binary systems, LMXBs. (LMXBs consist of a neutron star or black 
hole accreting from a companion less than one solar mass. 2 ) The recent detection of coherent X- 
ray pulsations with a millisecond period from a suspected low-mass X-ray binary system appears to 
confirm this link. 3 Here we report observations showing that the orbital period of this binary system 
is two hours, which establishes it as an LMXB. We also find an apparent modulation of the X-ray flux 
at the orbital period (at the two per cent level), with a broad minimum when the pulsar is behind this 
low-mass companion star. This system seems closely related to the "black widow" millisecond radio 
pulsars, which are evaporating their companions through irradiation. 4 8 It may appear as an eclipsing 
radio pulsar during periods of X-ray quiescence. 



The transient X-ray source SAX J1808. 4-3658 was 
first observed in September 1996 by the Wide Field Cam- 
eras on the BeppoSAX X-ray satellite during a bright 
state which lasted about 20 days. 9 The source was not 
detected (X-ray flux < 3 x 10 -11 erg cm -2 s -1 in the 
2-10 keV band) during an August 1996 observation, but 
reached a peak X-ray flux of 2 x 10 -9 erg cm~ 2 s _1 dur- 
ing the September bright state. Also during the bright 
state, BeppoSAX detected two type I X-ray bursts from 
the source, each lasting less than 30 s. Such bursts are 
due to thermonuclear ignition of accreted material on the 
surface of neutron stars that have low (< 10 10 G) mag- 
netic fields. 10 Analysis of the bursts in SAX J1808.4-3658 
indicates that it is 4 kpc distant and has a peak X-ray lu- 
minosity of 6 x 10 36 erg s _1 in its bright state and < 10 35 
erg s _1 in quiescence. 9 The source position was deter- 
mined to within a la radius of 0.6 arcmin. 

More recently, a serendipitous slew of the Proportional 
Counter Array (PCA) on the Rossi X-Ray Timing Ex- 
plorer (RXTE) over this region of the sky on 9 April 
1998 indicated the presence of an X-ray source, desig- 
nated XTE J1808-369. Repeated scans over the region 
with the PCA on 11 April localized the position with suf- 
ficient precision to confirm that the source was probably 
the same as SAX J1808.4-3658. 11 RXTE subsequently 
made numerous pointed observations of the source during 
11 April-6 May on a public "target-of-opportunity" basis. 
The source reached a peak luminosity of 5 x 10 36 erg s _1 
on 13 April and had faded below 10 35 erg s _1 by 6 May. 



Strong coherent 400 Hz pulsations were clearly de- 
tected in the 2-30 keV PCA data from most of these obser- 
vations, and were first reported by Wijnands and van der 
Klis. 12,3 These are the fastest persistent coherent pulsa- 
tions ever detected from an X-ray binary. From standard 
magnetic disk accretion theory, the detection of X-ray pul- 
sations at a luminosity ~ 10 36 erg s^ 1 requires that the 
neutron star's surface dipole magnetic field be B < 10 8 G 
(ref. 3). This is much weaker than the ~ 10 12 G fields 
found in other (slower) accretion-powered pulsars, 13 but 
it is consistent with the weak fields expected for type I X- 
ray bursters. 10 If pulsations are detected when the source 
is at lower luminosities, the implied surface field would be 
even weaker. 

In order to search for orbital Doppler shifts, we first 
selected the 3-30 keV photon arrival times during 11-18 
April and binned them into 2 -11 s (w0.5 ms) samples, af- 
ter correcting these times for RXTE's motion with respect 
to the solar system barycenter using the best BeppoSAX 
position for the source 9 (right ascension 18h 08m 29s, dec- 
lination —36° 58.6', equinox J2000.0). We then measured 
the barycentric pulse frequency at 128 s intervals using 
an oversampled Fourier power spectrum. The sequence of 
pulse frequencies showed an obvious 2 hr sinusoidal mod- 
ulation which was well fit by a constant spin frequency 
plus a circular, Keplerian orbit. 14 

Using this provisional timing model, we epoch-folded 
each 128 s interval in the 11-18 April data and cross- 
correlated the resulting pulse profiles with a sinusoid to 
measure the pulse phase history, which generally yields 
more precise timing parameters than a pulse frequency 
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Table 1. Observed parameters of SAX J1808.4-3658 
Barycentric pulse frequency, u (Hz) 400.9752106(8) 
Pulse frequency derivative, |z>| (Hz s _1 ) < 7 x 10 -13 

Projected semimajor axis, aisini (lt-ms) 62.809(1) 
Orbital period, P orb (s) 7249.119(1) 
Epoch of 90° mean longitude, T„ /2 (MJD) 50914.899440(1) 
Eccentricity, e < 5 x 10~ 4 

Pulsar mass function, /i (M ) 3.7789(2) x 10" 5 

Modified Julian date MJD = JD - 2400000.5. Numbers in 
parentheses are the la uncertainties in the last significant fig- 
ure, and upper limits are quoted at the 2a level. The v limit 
refers to the magnitude, independent of sign. Epochs are given 
in Barycentric Dynamical Time (TDB). 



history. These pulse phases were well fit (reduced x 2 =l-26 
with 441 degrees of freedom) by a constant spin frequency 
v plus a circular Keplerian orbit. The best-fit parameters 
are given in Table 1. The pulse arrival time delays with 
respect to a constant frequency model are plotted in Fig- 
ure 1. No improvement in the fit to the 11-18 April data 
is obtained by including an orbital eccentricity e or a spin 
frequency derivative v. A more sensitive search for a non- 
zero value of v, which includes the later data taken as the 
X-ray outburst faded from view, will be presented else- 
where. 

The present limits on v constrain the positional uncer- 
tainty of the star. An error of 0.6 arcmin in the assumed 
source position would give rise to periodic timing residu- 
als with 1 yr period and ^80 ms amplitude, 15 and such 
residuals would lead to an apparent v in the timing so- 
lution for observations shorter than 3 months. Our limit 
on v from the 11-18 April data suggests that the Bep- 
poSAX position is accurate to within 25 arcsec in ecliptic 
longitude. We note that a probable optical counterpart 
has been identified, lying 19 arcsec from the BeppoSAX 
position. 16 ' 17 

The pulsar mass function /i, which relates the pulsar 
mass mi, the companion mass ni2, and the binary inclina- 
tion i (the angle between the line of sight and the orbital 
angular momentum vector), may be computed from the 
observed Keplerian parameters a\ sin i and P orb , 
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from which we find f 1 = 3.7789(2) x 1CT 5 M Q (here Mq is 
the solar mass). Given assumed values of mj and i, equa- 
tion (1) can be solved for m 2 . For m 1; dynamical mass 
measurements in both binary millisecond radio pulsars 18 
and accretion-powered pulsars 19 are all generally consis- 
tent with a neutron star mass of 1.35 Mq. However, neu- 
tron stars in LMXBs have undergone substantial accretion 
from their binary companion and may be as massive as 2 
Mq (e.g., ref. 20). For i, we note that for an ensemble of 
binaries whose inclinations are distributed randomly, the 
a priori probability of observing a system with inclination 
i or smaller is (1 - cosz). If SAX J1808.4-3658 is drawn 
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Figure 1. Pulse arrival time delays due to the 2 hr binary 
orbit, with respect to a constant pulse frequency model. 
Mean longitude is the orbital phase angle measured from 
the ascending node, so that the pulsar is behind the com- 
panion at 90°. The squares are the measured time delays, 
the solid curve is the best-fit orbital model, and the dia- 
monds are the model residuals (on a 50 x expanded scale). 
The residuals have an r.m.s. uncertainty of 75 /js, shown 
for one of the points. 

from such an ensemble, then there is a 95% probability 
that m 2 < 0.14M Q if mi = 1.35 Mq. For mi = 2M , we 
find m 2 < 0.18 Mq. In either case, the binary separation 
is most likely of order 1 lt-s. 

The nature of the companion star provides an im- 
portant clue to the evolutionary history of an LMXB. 
The only direct information we have on the companion 
comes from observations of the probable optical counter- 
part, which suggest that the companion is a faint low- 
mass star subject to X-ray heating by the pulsar. 16 ' 17 
We can use the binary parameters to deduce more about 
this star. The companions in most LMXBs fill (or nearly 
fill) their critical gravitational potential surface, known 
as the Roche lobe. Because the mean density of a Roche- 
lobe-filling companion (for the case where m 2 < mi) is 
uniquely determined by the binary period, 21 we can relate 
the companion's radius in SAX J1808. 8-3658 to its mass 
as i?2 = 0.17 Rq (m 2 /0.1 Mq) 1 / 3 , giving a minimum com- 
panion radius of 0.12 Rq (here Rq is the solar radius). 
The assumption of Roche lobe overflow thus strongly con- 
strains the nature of the companion. Given the compact 
size of the binary, only white dwarf (WD) or low-mass 
main sequence (MS) companions are plausible. (In prin- 
ciple, a helium-burning star could also fit in this binary, 
but we would expect it to be much brighter than the pro- 
posed optical counterpart. 22 ) A WD with mass m has ra- 
dius R = 0.013 i? Q (l+X) 5 / 3 (m/M Q )- 1 / 3 , where AT is the 
hydrogen mass fraction. 23 Comparing this to our equation 
for i? 2 , we see that a helium WD (X = 0) would be too 
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small, and hence is ruled out. Even a hydrogen-rich WD 
(X = 0.9) is essentially excluded for all but i « 90°. How- 
ever, the lack of a deep X-ray eclipse rules out i > 82°, 
assuming a Roche-lobe-filling companion. 

Normal low-mass hydrogen MS stars have R/Rq « 
ui/Mq, and comparison of detailed models 24 with our 
equation for i?2 yields a Roche- lobe-filling mass of 0.17 
Mq. This would require a small binary inclination (i < 
20°), which has low a priori probability (5%). However, 
normal MS models are inappropriate in our case, since 
the stellar structure must be strongly influenced by irra- 
diation from the pulsar. If the incident flux at the com- 
panion surface exceeds ~ 10 10 erg cm -2 s _1 , the star may 
be "bloated" to larger radii. 25 This is especially true for 
MS stars having m < 0.3 Mq, which have fully convec- 
tive envelopes. For the ~1 lt-s binary separation in this 
system, the X-ray flux at the companion will exceed 10 14 
erg cm~ 2 s _1 , and the bloating can be severe. Under these 
conditions, MS stars as light as 0.1 Mq or less may fill their 
Roche lobe, consistent with a more likely inclination. 

Another possibility is that irradiation by the neutron 
star might directly drive mass loss from the companion 
("ablation") whether or not it fills its Roche lobe. 26 The 
binary parameters of SAX J1808. 4-3658 are very similar 
to those of the five "black widow" millisecond radio pul- 
sars in very close binaries, all of which are ablating their 
low mass companions. 4-8 At least four of these radio sys- 
tems show eclipses which are too broad to be caused by a 
Roche-lobe-filling companion, and which are instead at- 
tributed to an ablated wind. 4-6,8 

We find evidence that the X-ray flux from SAX 
J1808. 4-3658 is slightly modulated at the orbital period. 
The apparent modulation is roughly sinusoidal with 2% 
amplitude and a minimum when the pulsar is behind the 
companion, and it is much broader and much shallower 
than the (at most) 5 minute eclipse possible from a Roche- 
lobe-filling companion. However, the detailed features of 
this 2 hr modulation must be viewed with caution, since 
its strength is comparable to the variation of the RXTE 
background with the 96 min spacecraft orbit. The simi- 
larity in strength and time scales of these two variations 
makes them difficult to disentangle, and final confirma- 
tion must await either additional data or an improvement 
in the RXTE background model. Still, the detection of a 
2 hr source flux modulation (if not its precise morphology) 
seems secure. 

We suggest that the intensity dips are due to scatter- 
ing in an ablated wind. Given the similarity between this 
source and the eclipsing radio pulsars, SAX J1808. 4-3658 
may emerge as a radio pulsar during X-ray quiescence. 27 
Moreover, since the radio emission would be less pene- 
trating than the X-rays, a deeper eclipse dip might be 
observed, possibly providing a strong constraint on the 
binary inclination. The presence of the slight X-ray dips, 
if confirmed, would rule out small binary inclinations and 



make it highly likely that the companion mass is less than 

0. 1 Mq. 

Whether the mass accretion is fed by Roche-lobe over- 
flow or an ablated wind or both, we can understand the 
transient nature of the X-ray emission as long as some 
of this material forms an accretion disk. Although the in- 
stantaneous mass accretion rate during the 1996 and 1998 
bright states was « 3 x 10~ 10 Af© yr -1 , the small duty 
cycle of the X-ray activity indicates a long-term mean 
mass transfer rate of M pa 1 x 10~ n Mq yr -1 . For such 
a low M, the accretion disk would probably be subject 
to dwarf-nova-type instabilities, leading to episodic out- 
bursts of X-ray emission. 28,29 We note that for a Roche- 
lobe-filling companion, angular momentum losses due to 
gravitational radiation would drive a mass transfer rate of 
10 1 Mq yr -1 for a 0.05 Mq secondary, 30 consistent 
with our inferred inclination constraints. 

Future X-ray outbursts from SAX J1808.4-3658 may 
allow detection of a spin frequency derivative, which would 
further constrain both the mass accretion rate and the pul- 
sar's magnetic field strength and would provide a probe 
of magnetic disk accretion torque theory in the previously 
unexplored regime of a very small magnetosphere. The 
characteristic accretion torque expected for the observed 
X-ray luminosity should cause a u of order 10~ 14 Hz s _1 , 
which is detectable in observations spanning a month or 
more. Further observations may also detect orbital pe- 
riod evolution, which would probe the competing effects 
of gravitational radiation, tidal interactions, and mass loss 
on the binary angular momentum. An astrophysically in- 
teresting limit of |P or b| < 10~ n could be reached in less 
than a year of observations. 
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